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Abstract 

Recently, there have been studies of parametric resonance decay of oscillating real homogeneous 
cosmological scalar fields, in both the narrow-band and broad-band case, primarily within the 
context of infiaton decay and (p)reheating. However, many realistic models of particle cosmology, 
such as supersymmetric ones, inherently involve complex scalar fields. In the oscillations of complex 
scalars, a relative phase between the oscillations in the real and imaginary components may prevent 
the violations of adiabaticity that have been argued to underly broad-band parametric resonance. In 
this paper, we give a treatment of parametric resonance for the decay of homogeneous complex scalar 
fields, analyzing properties of the resonance in the presence of out of phase oscillations of the real 
and imaginary components. For phase-invariant coupling of the driving parameter field to the decay 
field, and Mathieu type resonance, we give an explicit mapping from the complex resonance case to 
an equivalent real case with shifted resonance parameters. In addition, we consider the consequences 
of the complex field case as they apply to "instant preheating," the explosive decay of non-convex 
potentials, and resonance in an expanding FRW universe. Applications of our considerations to 
supersymmetric cosmological models will be presented elsewhere. 
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1 Introduction 



With the advent of inflationary theories of the early universe, it has been argued that the present 
stage of hot FRW "big-bang" cosmology was preceded by an epoch of cosmological evolution domi- 
nated by the dynamics of scalar fields The success of inflationary models in providing explana- 
tions for flat large-scale geometry (as suggested by the location of the acoustic peaks in the CMBR 
anisotropics), and for the origin of approximately scale-free adiabatic density perturbations (which 
can be used to simultaneously fit both the CMBR anisotropics and observations of cosmic structure 
formation), lends support to the idea of an early scalar- field dominated epoch. A crucial question 
in this picture is the nature of the transition from the scalar field dominated epoch, to the hot FRW 
epoch, which is referred to as reheating. The nature of this transition also relates to other aspects of 
early universe dynamics necessary for a successful cosmology, such as mechanisms of baryogenesis, 
the resolution of cosmological moduli problems, and possible sources for non-thermal dark matter. 

The standard approach to reheating, which applies to sufficiently weakly coupled inflaton fields 
0, is to treat quanta of the infiaton field as particles, which undergo independent single particle 
decay; this treatment, if adequate, has the advantage that the post-infiation reheat temperature is 
determined by the microphysics of the model. For infiaton fields with mass as suggested by the 
simplest chaotic or supersymmetric models, and decaying by gravitational strength interactions, 
this treatment is adequate, leading to moderate reheat temperatures (< O(10^°) GeV) which are 
consistent with the absence of GUT-scale defects such as monopoles, which are capable of incorpo- 
rating a variety of (s)leptogenesis or electroweak mechanisms for generation of the BAU, and which 
avoid, in the supersymmetric case, cosmological problems with thermal overproduction of gravitinos 
after reheating 0. 

Recently it has been realized that the standard treatment of reheat in terms of single particle 
infiaton decay may be seriously misleading in circumstances where there is coherent enhancement 
of the transition to bosonic decay products 0-0. For large mode occupation numbers of the decay 
product field we may treat its dynamics as being essentially classical. Mode by mode for the decay 
product field its coupling to the oscillating infiaton field induces a periodic time dependence in the 
mode mass (modulated by cosmic expansion which "sweeps" the time dependence of each comoving 
mode through the bands of the stability chart of the mode equations.) This periodic modulation of 
the parameters of the oscillator associated with each mode of the decay product field, can induce 
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parametric resonance in bands of the mode parameters, leading to exponential growth in the decay 
mode amplitude. The resonant decay of the inflaton may have important cosmological implications 
like non-thermal symmetry restoration and subsequent formation of topological defects [T^,|T^, 
revival of GUT baryogenesis scenarios |TP,2D|, supersymmetry breaking |]2T]], superheavy particle 
production ||T^,|22|, and gravitino production 



The exponential growth in the mode occupation number may be modified or regulated by a 
number of physical processes. These include the decay of produced quanta to other particles 
^ or the rescattering of final state particles [Q. Another possibility, occurring in models with 
gauge-strength self-interactions between the produced final state particles, is the regulation of the 
parametric resonance by the self-interaction induced effective masses of the produced quanta, which 
can move these quanta out of the available resonance bands; in this case, resonance only proceeds 
as thermalization and Hubble dilution reduce the plasma masses of the final state quanta, resulting 
in a quasi-steady-state resonance conversion of inflaton oscillation energy to decay products ||2^ . 
This general scenario for the physically realistic case of decay products with gauge charge has been 
verified in explicit calculations in the narrow-band resonance case [^, and in numerical simulations 
in the broad-band case ||10|| . 

While to date analytical and numerical treatments of parametric resonance have considered 
oscillations of a single real field decaying to a single real decay product field, in realistic models the 
field content is often more extensive. In the case of supersymmetric theories this is unavoidably 
the the physical scalars of simple (N=l) supersymmetry come as components of chiral 

multiplets and are complex. So for these types of theories, we should at the very least consider the 
nature of coherent decays when the fields involved are complex, though non-supersymmetric models 
with multiple real scalar fields may share some of the features of the simplest complex case. 

Within supersymmetric models of particle physics, there are several different circumstances 
under which the decay of a homogeneous complex scalar condensate may occur in the early universe. 
At the end of inflation one expects to have a spatially homogeneous inflaton scalar condensate, whose 
decay energy will ultimately be responsible for cosmic reheating and the initiation of hot big-bang 
cosmology. As well, in the supersymmetric standard model there are directions in the scalar fleld 
space of squarks and sleptons which are F-flat and D-flat, and which only gain a potential from 
supersymmetry breaking. In the early universe these directions may be populated with (very) large 
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vev's after the end of the inflationary epoch; these vev's may carry enormous vev to mass ratios 
(Mathieu resonance parameter q) and couple to other directions in scalar field space with couplings 
capable of inducing resonant decay. Finally, supersymmetric models are generically plagued with 
gauge singlet scalar moduli, whose homogeneous oscillation poses grave cosmological difficulties 
which might be ameliorated by coherent decay of their oscillation amplitude. 

For self-interactions of complex scalars of the general form dictated by the F-term and D- 
term couplings arising in globally supersymmetric theories, the fields generally appear in complex 
conjugate pairs for the F-terms and diagonal D-terms. For example, let us consider a complex 
scalar field $ in a chiral supermultiplet whose decay will be induced by a trilinear (renormalizable) 
couphng in a superpotential 1^ to a chiral multiplet labelled by its scalar S, where W ~ $SS. The 
resulting F-term coupling inducing the decay is then of the form $*$S*S, and is invariant under 
global phase redefinitions of either the $ or the S. We will see in the next section that in cases 
such as this the phase invariance of the resulting couplings implies that the equations for modes 
of the real and imaginary components of S are decoupled and independent, and will allow us to 
simply analyze the resonant decay of a $ condensate with out of phase oscillations for the real and 
imaginary components of into real and imaginary components of the decay product field S. 

We can always phase rotate our scalar field $ to a basis such that its initial vev lies along the 
real axis. If there is no component of force along the direction of the imaginary axis (i.e. the scalar 
potential is phase invariant), the trajectory of the motion of $ is limited to the real axis and the 
field hits the origin as it oscillates back and forth. In this case, provided that the coupling of the 
oscillating field to the final state field is also phase invariant, the situation is exactly that of a 
real oscillating field, and the same arguments apply for parametric resonance particle production. 
However, if the scalar potential is not phase invariant, i.e. depends on the phase of the oscillating 
field as well, a torque is exerted on the field. This leads to the defiection of the trajectory from a 
straight line and results in changing the trajectory into something that finally resembles an ellipsoid, 
after the torque in field space has effectively ceased its action. In this case the field no longer passes 
through the origin but rather has a finite distance of closest approach to it. This will have important 
implications for broad-band parametric resonance as we discuss below. 

In supersymmetric models not only are complex scalar fields inherently involved, but also a 
phase dependent part of the scalar potential can arise naturally from supersymmetry breaking. 
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Let us consider the simplest case with the following terms only involving the inflaton in the su- 
perpotential W = |m$^ + |A$'^. In supergravity models with broken supersymmetry, there is a 
corresponding phase dependent term (the "A-term") Am^^-+h.c. in the scalar potential, where A 
is a dimensionless model-dependent constant and mA is the scale of supersymmetry breaking in the 
sector in which $ lies. There is also a phase dependent term mA*$$*^ + h.c. in the F-term part of 
the scalar potential. This generically occurs in minimal supergravity models for inflation where the 
superpotential contains a series of \nj^^ terms and occurs even in "no scale" supergravity 
models after the inclusion of radiative coreections to the effective potential p8 |. 



For V D the potential along the angular direction is periodic with m — n minima. In 

general, during inflation $ rolls down to its minimum both along the radial and angular directions. 
In order to have a torque to deflect the trajectory, $ must not be at the minimum along the angular 
direction at the onset of radial motion. This can happen in two ways: either there are several 
phase dependent parts of the potential with a non-adiabatic transition from the minimum of one 
to another, or $ does not settle at the minimum along the angular direction. The first possibility 
happens when other supersymmetry breaking sources in the early universe (e.g. non-zero energy 
density of the universe or finite temperature effects) are dominant over the low energy one. In this 
case the minimum in the angular direction at early times is different from that at late times (due to 
independent phases for the coefficients of different A-terms). If the transition from one minimum to 
another one is non-adiabatic, $ will not be at the minimum at late times regardless of its start at 
the minimum at early times. The second possibility happens when the potential along the angular 
direction is fiat enough during inflation. In this case $ will not roll down to its minimum and can 
start at any position at the onset of radial motion. In both cases, the further $ is away from the 
minimum, the larger the deflection of its trajectory and the wider the ellipsoidal shape will be. 



2 Complex Mathieu Resonance 

As described in the previous section, the potential for complex scalar oscillations usually includes (as 
well as the scalar mass-squared terms) Hubble induced A-terms which are mainly effective during 
the first few cycles of oscillation. The A-terms provide a "torque" to the complex oscillation during 
the first few cycles, resulting in a net "elliptic" motion in the mass-term induced potential, after 
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the A-terms have effectively ceased to be active. The resulting elliptic oscillation in the mass term 
potential will be damped by the Hubble drag, resulting in the ellipse shrinking over time. 

In order to introduce new considerations characteristic of resonance with complex fields, without 
getting involved in model dependent details, in most of this paper we shall simply ignore the damping 
and consider complex, elliptic, constant amplitude oscillations. In particular, this means we do not 
need to specify which particular field is considered (e.g. infiaton versus susy standard model fiat 
direction), nor do we need to determine the cosmological details involved in determining expansion 
and damping at the time that the oscillations of the field in question occur. In addition to presenting 
a tractable and interesting mathematical problem, consideration of the undamped oscillation should 
also provide the essential features of the full cosmological case including the effects of expansion, 
at least in the generic case of broad-band resonance. This follows from the key observation of 
that in the broad-band case the resonant excitation of the decay product field occurs over a tiny 
fraction of the cycle of the driving field, when the latter passes near the origin, as only here do 
the decay product field dynamics depart from the adiabatic regime. So mode number excitation 
proceeds by a series of abrupt jumps, and the dynamics of a given jump may be considered with the 
instantaneous value of the oscillator parameters, resulting in the picture of "stochastic resonance" 
analyzed in [Q. The present paper discusses the changes in the dynamics of decay mode excitation 
which arise from the complex nature of the driving field oscillation — the differences in question arise 
from suppression of the adiabaticity violations that induce the jump in mode occupation numbers 
of the decay product field, so we expect that considerations using the instantaneous value of the 
driving oscillation amplitude should give insight in the complex case, much as such considerations 
did in the real stochastic resonance case. 

In any case, for the purposes of our present calculations we shall consider the parametric res- 
onance production of decay product field modes H, from phase-invariant coupling to constant- 
amplitude out of phase ("elliptic") oscillations of a driving field $. Detailed cosmological studies 
of applications to infiaton or moduli oscillations will be considered elsewhere. 

With the couplings discussed in the previous section, after the A-terms cease to be effective, the 
equation of motion for the S field is of the form: 
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where is the decay product field mode with comoving wavenumber k, a is the FRW scale factor, 
m= the mechanical mass of the S, and the superpotential coupling is as above. We note that both 
the real and imaginary piece of the S field will separately obey this equation, and hereafter we use 
X to denote either the real or imaginary piece of S. 

In our analysis, we will treat the physical momentum of the decay field mode and the relative 
phase and amplitude of the driving field oscillation as fixed parameters, and attempt to map out 
the regions of instability in their parameter ranges. As noted above, for the case of stochastic 
broad-band resonance, where the amplification occurs in small intervals while the field passes close 
to the origin, one should be able to approximate the instantaneous behaviour within each interval 
by the corresponding behaviour of a system of the type we analyze here. 

We decompose the driving field $ into real and imaginary pieces as follows: 

$ = + ^0/- (2) 

By a phase rotation we put the largest amplitude component of oscillation into the real piece, and 
so we write: 

(j)R = 0sin(m^t) (3) 
0/ = f(j)cos{m^t), (4) 

where now (j) is the constant amplitude of the real component of oscillation and / G [0, 1] is the 
"out of phase" fractional component giving elliptic oscillation in the complex $ plane; we will be 
particularly interested in the case where / <^ 1. 

We wish to cast this into the canonical form of the (real) Mathieu equation: 

y" + {A-2qcos{2z))y = 0, (5) 

where ' denotes derivative with respect to the independent variable z. We begin by substituting 
our definition of the $ field into (|ID above, giving 



Xk + {^ + ml + g^4>'^ (sin2(m^t) + f cos^{m^t))^ Xk = 0. (6) 



We replace cos (m^pt) with 1 — sin (m^t) and collect terms, giving 

Xk+(^ + ml + fg^' + (1 - f) g^' sm\m,t) ) Xk = 0- (7) 
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Using the half-angle formula sin 




(1 — cos 2^) we obtain the form 



Xk+('^ + ml + fg^4>' + - /') 9'4>' (1 - cos{2m^t))) Xk = 0. 



which may be rewritten in the form of the Mathieu equation 



xl + iAkif)-2qif) cos22;)xA: = 



(9) 



with the following new identifications 



z 



(10) 



?(/) 



ml 

(1 - fW 



+ Hf) 



(11) 



(12) 



Ami 



Notice that the coefficients of the Mathieu equation are now functions of the imaginary fraction 
/. This is an important feature, as it means that the characteristic behaviour of the parametric 
resonance as described by the Mathieu equation takes the same form in the complex case as in the 
real case; however, the relationship between the physical parameters of the process and the Mathieu 
coefficients is redefined. 

So we see that there is a mapping that takes the Mathieu equation for the complex modes of 
the decay field when driven by the complex parametric field with out of phase real and imaginary 
pieces in its oscillation amplitude, and maps it to a real Mathieu equation for the oscillations of 
the real and imaginary pieces of the decay product field with shifted parameters. There are several 
features of this mapping that simply encapsulate physical features of the original problem. 

First we note that for the case where / = 0, where the oscillation of the parametric driving 
field is along the real axis, the coefficients and q reduce to those previously considered in the 
literature for the case of purely real parametric oscillation [§, In the other extreme limit, when 
/ = 1, we have g(l) = 0, meaning that one is restricted to be along the Ak axis on the Mathieu 
equation stability diagram, allowing only non-resonant particle production. This corresponds to the 
physical observation that because the decay couplings were phase invariant, our original equation 
for the complex oscillations involved only the magnitude of $ in the oscillation equation. In the 
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case that the real and imaginary pieces of the $ oscillation have the same amplitude (/ = 1), 
then the coefficients in the x equation of motion become time independent and there can be no 
parametric amplification. We also note that the imaginary fraction / of the oscillations enters 
into the coefficients only as meaning that the effect is the same regardless of direction traveled 
around the ellipse (i.e. whether / is positive or negative). This reflects that fact that the original 
equation for the complex oscillations is second order and symmetric under time inversion, which 
corresponds to having the parametric fleld circulate about the oval in the $ plane in the opposite 
sense, or reversing the sign of /. 

Finally, in the intermediate regime / e (0,1), q{f) is always lessened, while is always 
increased (compared to / = 0), so increasing / never causes the system to leave the physical 
regime. As noted above, increasing / means moving "inland" on the stability diagram for the 
Mathieu equation. In general, this causes suppression of the resonant growth of the x modes; 
however, it also allows one to explain the counterintuitive observation that in certain cases the 
resonant band exponential growth parameter /ik may actually increase as one turns on the out of 
phase component /. To understand this, imagine that for oscillations with no imaginary fraction 
/ one is sitting in parameter space at the lower border of one of the instability bands (where /ik is 
zero). Now slowly increase /. The parameter mapping derived above implies that you start to move 
in a "northwest" direction on the band chart into the instability band. As such, your /ik begins 
to increase. As you continue to increase /, you will eventually hit the maximum possible /ik along 
your trajectory, after which your /ik begins to drop again. Eventually you will leave the instability 
band altogether for sufficiently large /. For very high order instability bands, it should be possible 
to encounter many bands along one trajectory of increasing /. 

From a different point of view, were one to look at the instability diagram as a function of the 
the Ak and q of standard real Mathieu parametric resonance (ie. as a function of our Afe(O) and 
q{0), for different values of /), the effect of turning on / would be seen to manifest itself as both a 
narrowing and a downwards shift of the instability bands as / increased. In addition, isocontours 
of Hk would be seen to "flow out" of the bands as / increases. For f — 1 each instability band 
collapses to a line with /ik — 0, as with a phase independent coupling of the parameter fleld to the 
mode fleld there would be no time dependence in the mode fleld equation of motion, and the system 
would be stable for all Ak and q. Figure 1 illustrates in detail both the bending of the bands and the 
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decrease of /Xfc for the first resonance band as we turn on the imaginary fraction / of our oscillations. 
We expect that the suppression of resonance for fixed non-zero imaginary fraction / is stronger in 
higher resonance bands, as at larger q the resonance proceeds by violation of adiabaticity in the 
decay mode evolution, and for large q at fixed / the induced decay mode field mass is always large. 
This is illustrated in Figure 2 where we show the quenching of resonance for the first 7 resonance 
bands, as / is turned on. We see that there is more severe suppression for the higher bands, in 
accord with our physical intuition; in the next section we will analytically estimate the extent of 
the domain of surviving resonant bands, for fixed imaginary fraction /. 



3 Parameter Domain For Broad-Band Resonance 

Let us first briefly discuss the effect of mixing in the narrow-band regime. The narrow-band reso- 

4 

nance in the case of a real oscillating field is efficient for (y) ^ < < ^ [g]. In the case of a complex 
field with mixing parameter / this reads as 
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^ ^ ^ ^ < </> < ; , (13) 



(1 - P)g^ ) - - V (1 - P)g 

where m and g are replaced with ^Jvrfi -|- f'^g'^<\P' and a/1 — respectively. It is easily seen that 
the condition for an efficient narrow-band resonance remains almost unchanged in the complex case, 
unless / is close to 1. Therefore, in physically interesting situations the narrow-band resonance will 
not be substantially affected by the mixing. Of course, in the extreme case with / = 1, there is no 
time variation in the Mathieu equation and, hence, no narrow-band resonance. 

In the case of real broad-band parametric resonance, Kofman, Linde, and Starobinski 0, ^ argue 
that the requirement of adiabaticity violation for broad-band parametric amplification implies that 
it only occurs with significant /i^ for — 2g < y^g. Their argument presupposes that the decay 
terminates after it has entered an "explosive" phase, where the effective mass of the decaying is 
dominated by the coupling to the plasma of decay product x modes which has been built up by 
parametric resonance decay. The effective physical 3-momentum of the quasi-relativistic decay x 
modes is of order their energy, which is no more than of order the induced mass of the decaying 
0; this is of order (^Xend which in turn is of order (yf^end which can also be written y^^m^0^. So 
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(^phys/''^^) ~ which caii be rewritten as — 2q < ^g. For a detailed discussion we 

refer to the treatment in 

We have seen in the preceding section that the case of complex oscillation with imaginary fraction 
/ can be mapped onto a Mathieu equation with shifted resonance parameters. By substituting the 
"shifted" parameters induced by the non-zero imaginary component of oscillation, we should thus 
be able to determine what range of q supports broad-band resonance for oscillation with a given 
fraction of out of phase imaginary component for the oscillation of the driving parameter. 

Recall the expressions for the equivalent shifted Ak{f) and q{f) from equations ( pT]) and ([T^) 
respectively. Substituting these expressions into the relation — 2q < ^ allows us to write it as: 



(14) 



This leads us to the relation 

A,{0) - 2g(0) + Afq{0) < v^W^v^, (15) 

or, defining Ek = Ak{0) — 2g(0), we have 

Ek + 4f g(0) < ./l^y^). (16) 

If we recall that physical values of are positive semi-definite, we find that for a fixed non-zero 
imaginary fraction / there is an upper bound on the parameter g(0) for which resonance occurs, 
and the allowed range of resonant g(0) is bounded above as ^^jr- (For the small imaginary fractions 
/ of physical interest, however, the weaker approximate bound of j^f'^ will suffice). So instead 
of an an ever widening resonance region above the Ak = 2q line, with thickness of order ^/g, as 
one has in the real case, in the complex case with fixed non-zero imaginary fraction /, one instead 
has a region above the Aj^ = 2q line of finite extent, with an upper bound on the values of the 
q parameter which can result in resonance. This is qualitatively reasonable, as a fixed imaginary 
fraction / for the oscillation means that as we scale up q the ellipse of $ broadens as it lengthens, 
preserving its shape; so throughout the $ oscillation |<l>| has a large value, inducing a large mass for 
the modes of the decay field S, which in turn leads to adiabatic evolution of the H, and suppression 
of broad-band parametric resonance production of the S. 
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4 Complex Resonance in "Instant Preheat" 



Recently, a simpler method of efficient scalar field decay has been proposed, called "instant preheat" 
[]16| . Within models of this type the decaying field rolls once through the origin, at which point 
the mass of the decay product field to which it is coupled passes through zero, and modes of the 
decay product field experience non-adiabatic excitation. As the decaying field rolls away from zero 
(perhaps monotonically) the modes of the decay product field grow in mass; they drain energy from 
the decaying field through their mass. As the mass of the modes of the decay product field grows, 
so does their decay width; their subsequent decay, after their mass and decay width have grown 
sufficiently, then releases the energy they have taken from the original decay field, and dumps it 
into their final decay products, which thermalize the resulting energy. 

It is interesting to note that while final state effects such as rescattering, backreaction, or plasma 
masses can prevent preheating from occurring, they are unimportant in the instant preheating 
scenario. The reason is that for these effects to become important, (at least) several oscillations 
are needed to build up a large enough occupation number for the final state field. In the instant 
preheating, on the other hand, the energy drain from the oscillating field occurs during each half 
of an oscillation. In fact, instant preheating can be efficient even if the adiabaticity condition is 
violated only during the ffist half of the ffist oscillation. Therefore, instant preheating is essentially 
unaffected by the final state effects. The mixing of the real and imaginary parts of the oscillating 
field, on the other hand, has the same effect in the instant preheating case as in the standard 
preheating scenario. We recall that the torque from A-terms deflects the trajectory of the oscillating 
field from that of a straight line. The Hubble induced A-terms have their largest value at the 
beginning of the oscillations, and rapidly decrease with Hubble expansion. This means that the 
deflection is largest during the initial oscillations. Thus, a large enough / to restore adiabaticity in 
the preheating case could do the same for the instant preheating case. 



5 Non- Convex Potentials 

Another possibility to achieve rapid decay of a homogeneous condensate occurs in the case that 
the potential governing the evolution of the condensate scalar has non-convex behaviour over some 
region of field space [0, 53]. In this circumstance, it becomes energetically favorable for a scalar 
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condensate in the non-convex region to decompose into inhomogeneous modes; provided the in- 
homogeneity occurs over long enough wavelengths, the price one pays in kinetic energy for the 
inhomogeneity is more than compensated by the decreased average potential energy of the regions 
of field excess and deficit compared to the average field value. This produces a wavenumber band for 
exponential growth of the mode amplitudes. This has been considered in both the case of inflaton 
decay |jl2|], and in the case of the growth of inhomogeneities in scalar condensates corresponding to 



F-flat and D-flat directions of the standard model with non-convex potentials of the type arising 



from gauge-mediated supersymmetry breaking ||29|. It is clear that this is one type of instability 



which is not vitiated by having the scalar order parameter complex or involving multiple scalar 
fields. If there is a region in field space with respect to which the potential is non-convex in some 
direction, then fluctuations corresponding to modes of the field variation in that field direction, 
of sufficiently long wavelength, will win on the potential versus kinetic energy budget, and grow 
exponentially. Indeed the treatment of (complex) flat directions in the supersymmetric standard 
model in [^] explicitly analyzes the conditions for instability of a complex field with a potential 



which is a non-convex function of the field modulus, and exhibits the resulting instability bands. 

6 Other Couplings 

Here, we briefiy comment on the situation for another type of coupling between $ and S fields 
which is also of interest and application. This is the coupling g'^{4>RXR + 4>iXi)^ where its simplest 
manifestation is for the potential ^($) = ;^A|$|'^, with $ and S being the same field. It also arises 
in supersymmetric models from the D-term part of the scalar potential. This type of coupling leads 
to the mixing of Xr and Xi mode equations: 

XR,k + + m^^g'^<pR^^ XR,k + 9'^<pR(t>iXi,k = (17) 
Xi,k + + ni^^g'^(j)i'^^ xi,k + g'^(t)R4>iXR,k = 0. (18) 

In this case the mass eigenstates are and instead of (pn and 0/ themselves. For 

oscillatory motion of (pR and (pj with a phase difference, there are two periodic changes that may 
lead to resonance: change in the mass eigenvalues (the usual parametric resonance) and change in 
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the mass eigenstates. They can't be simply superimposed and it is not very easy to give rough 
arguments for the instabihty bands and the respective value of fikS. The important point is that 
for such a coupling, even in the / = 1 case there is still time variation in mode equations. This 
variation is present in both the mass eigenstates and mass eigenvalues. 

7 Cosmic Expansion and Complex Resonance 

So far, we have considered modifications to parametric resonance decay which arise in complex field 
oscillations in the absence of effects of Hubble expansion. As we noted above, since broad-band 
resonance is induced by non-adiabaticity of the x evolution during small intervals of the oscillation, 
instantaneous approximation of the x excitation should be a useful guide during each of the jumps 
in mode number. Cosmic expansion then functions to shift the parameters of the oscillator between 
episodes of mode excitation as (p passes near zero. We now examine the implications of this in both 
the narrow- and broad-band cases. 

Implications for the narrow-band case are simple; as we have seen, the introduction of a phase 
difference between real and imaginary components of our complex inflaton field $ only kills narrow- 
band resonance for phase differences approaching |, or, in the language of this paper, for / = 1. 
Therefore, the resonance should be qualitatively the same in the static approximation and with the 
Hubble expansion included. 

For broad-band resonance the situation is completely different. According to equation ([16D , the 
broad-band resonance is shut off for q > jqI'^- In the static limit / and q are both constant 
and resonance is either suppressed, or viable. In an expanding universe, / eventually becomes 
approximately constant as the Hubble induced A-terms turn off (indeed, as pointed out earlier, 
after several Hubble times the motions along the real and imaginary directions are decoupled and 
free), while, on the other hand, q{t) = (^^) is redshifted as a~^. This implies that even if the 
resonance is suppressed initially, it may be initiated after a sufficient time such that q(t) < j^f'^- 
The right-hand side is less than 1 (or very close to it) for / ^ |. Therefore, in the case of large out 
of phase components of oscillation, the broad-band resonance is killed and resonance may resume 
only in the narrow-band regime at a later time. In most physically interesting cases, however, 
f < \ and the right-hand side is considerably greater than 1. In such cases, broad-band resonance 
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is not eliminated in an expanding universe, but rather delayed. The parameter / is determined by 
the action of the scalar potential (including A-terms) for the oscillating field, and after the initial 
oscillations it often becomes effectively time-independent. Depending on the dimensionality of the 
A-term, it may be a function of g^, the value of q at the start of oscillations. If < ]^/~^, the 
onset of broad-band resonance will be unaffected. For > -^f'^, resonance is delayed initially, 
but will resume after sufficient expansion such that q < qeq = A larger / leads to a smaller 

qeq, for / ^ I we have geq < 1 and resonance can only occur in the narrow-band regime. For / = 1 
resonance is eliminated. 

Even though the broad-band resonance (for interesting cases) is only delayed in an expanding 
universe, the mixing still has important consequences. Perhaps the most notable example relates 
to the production of superheavy particles during resonance. In the standard preheating, S's with a 
mass up to q^m^ can be produced. A reduction in q at the onset of resonance implies a reduction in 
the maximum mass of produced particles. This is even more pronounced in the instant preheating 
case. Here S decay products \1/ with masses up to q^m^ and which have a large enough coupling 
h to S, can be produced ||16[. A smaller geq has a two- fold effect in this case. First, the allowable 
masses are smaller, and second, the decay rate F^ = ^g(p may not be large enough (compared to 
the frequency of oscillations m^) for efficient production of \I''s. It is easily seen that F^^ < for 
h < Att^ f. Therefore, \l/ production is not efficient h < Att^ f. Even for h ^ An^ f, only \I''s with 
a mass < -^m^ can be produced. 

8 Conclusions 

In this paper, we have considered the changes to the standard picture of parametric resonance decay 
of a real homogeneous cosmological scalar field which arise if the field is instead complex, with out 
of phase oscillation of its real and imaginary components and a phase invariant decay coupling. For 
the case of complex Mathieu type resonance, we give an explicit mapping to a corresponding real 
Mathieu resonance with shifted parameters that encode the effects of the out of phase components of 
the oscillating decay field. We showed the resulting effects on the instability bands, demonstrating 
how they shift and shrink with increasing out of phase ("elliptic") component of the driving field 
motion, limiting the swath of instability to a finite area on the A^-q chart, and eliminating broad- 
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band resonance in the higher modes. We argued that similar effects may be present in the case of 
complex field models of "instant preheat," but that instabilities due to regions in field space with 
non-convex potentials are qualitatively the same in the complex case. Finally, in the context of an 
expanding FRW universe, we noted that the presence of a fraction of out of phase oscillation would 
usually delay the onset of parametric resonance, but not eliminate it entirely. 
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Figure Captions 



Figure 1: The Mathieu equation stability diagram for the first resonance band. The physical region 
lies above and to the left of the diagonal line A^. = 2q. Contour lines represent isocontours of n 
starting from fi = (band boundary) and increasing by units of 0.1 as one moves from left to right 
in the band. Different panels represent different imaginary fractions /: (a) / = 0; (b) / = 0.2; (c) 
/ = 0.4; (d) / = 0.6. 



Figure 2: The Mathieu equation stability diagram for the first 7 resonance bands. The physical 
region lies above and to the left of the diagonal line Ak = 2q. Contour lines represent isocontours 
of fi starting from fi = (band boundary) and increasing by units of 0.1 as one moves from left 
to right in the bands. Different panels represent different imaginary fractions /: (a) / = 0; (b) 
/ = 0.2; (c) / = 0.4; (d) / = 0.6. 
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